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The isoscalar and isovector particle densities and the surface tension coefficients at 
the average binding energy are used to derive analytical expressions of the neutron 
skin thickness and the isovector stiffness of sharp edged proton-neutron asymmetric 
nuclei. For most Skyrme forces these quantities are significantly larger than the well 
known ones. Using the analytical isovector surface energy constants in the framework 
of the hydrodynamical and the Fermi-liquid droplet models the mean energies and 
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the sum rules of the isovector giant dipole resonances are in fair agreement with the 
experimental data. 
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I. INTRODUCTION 

A simple and accurate solution of the particle density distributions was obtained within the 
nuclear effective surface (ES) approximation in Refs. [1-3]. It exploits the saturation properties 
of nuclear matter in the narrow diffuse-edge region in finite heavy nuclei. The ES is defined as the 
location of points with a maximum density gradient. An orthogonal coordinate system related 
locally to the ES is specified by the distance £ of a given point from this surface and tangent 
coordinates r\ parallel to the ES (see Fig. 1). Using nuclear energy density functional theory 
the variational condition derived from minimizing the nuclear energy at some fixed integrals of 
motion is simplified in the £, r\ coordinates. In particular, in the extended Thomas-Fermi (ETF) 
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approach [4], it can be done for any fixed deformation using the expansion in a small parameter 
a/R ~ A -1 / 3 <C 1 for heavy enough nuclei where a is of the order of the diffuse edge thickness 
of the nucleus, R is the mean curvature radius of the ES, and A is the number of nucleons. The 
accuracy of the ES approximation in the ETF approach without spin-orbit (SO) and asymmetry 
terms was checked [3] by comparing results with those of the Hartree-Fock (HF) and ETF theories 
for some Skyrme forces. The ES approach [3] was also extended by taking into account the SO 
and asymmetry effects [5]. 

In the present work, solutions for the isoscalar and isovector particle densities and energies in 
the ES approximation of the ETF approach are applied to analytical calculations of the neutron 
skin and isovector stiffness coefficient in the leading order of the parameter a/R (see also Ref. 
[6]). Our results are compared with older investigations [7-10] in the liquid droplet model (LDM) 
and with more recent works [11-19]. We suggest to study also the splitting of the isovector giant 
dipole resonances into the main and satellite (pygmy) peaks [18, 19] as a function of the analytical 
isovector surface energy constant of the ES approach within the Fermi-liquid droplet (FLD) model 
[20-22]. The analytical expressions for the surface symmetry energy constants are tested by the 
mean energies of the isovector giant dipole resonances (IVGDR) within the hydrodynamical (HD) 
and FLD models. 

The manuscript is organized as follows: In section II we give an outlook of the basic points of 
the ES approximation within the density functional theory and main results for the isoscalar and 
isovector particle densities. Section III is devoted to analytical derivations of the symmetry energy 
in terms of the surface energy coefficient, the neutron skin thickness and the isovector stiffness. 
The discussions of the results are given in section IV and summarized in section V. Some details 
of our calculations are presented in Appendixes A, B and C. 

II. ENERGY AND PARTICLE DENSITIES 

We start with the nuclear energy as a functional of the isoscalar and the isovector densities 

P± = Pn± Pp- 

E = J dr£( P+ ,p_) , (1) 
in the local density approach [4, 23-27] with the energy density £ (p + ,p_), 



£ (P+,P-) ~ ~b v P+ + Jj2 P+ + P+ i £ +(P+) -£-(/»+, P-)] + 
+ (C + + V +P+ ) (Vp + ) 2 + (C_ + V_ P+ ) (Vp_) 2 , 



(2) 
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where / = (N — Z)/A is the asymmetry parameter, N = f drp n (r) and Z — J drp p (r) are the 
neutron and proton numbers and A = N + Z. As usual, the energy density £ in Eq. (2) contains 
the volume part given by the first two terms of Eq. (2) and the surface part including the density 
gradients [1, 3]. The particle separation energy b v ~ 16 MeV and the symmetry energy constant 
of the nuclear matter J 30 MeV specify the volume terms in Eq. (2). Eq. (2) can be applied in a 
semiclassical approximation for realistic Skyrme forces [23-27], in particular by neglecting higher 
h corrections in the ETF kinetic energy [2-4] and Coulomb terms. Up to small Coulomb exchange 
terms they all can be easily taken into account (see Refs. [1, 3, 5]). The constants C± and T>± 
are defined by the parameters of the Skyrme forces [23, 24]. The the isoscalar part of the surface 
energy- density which does not depend explicitly on the density gradient terms, is determined by 
the function e + (p + ) [3, 5] which satisfies the saturation condition: e+(p) = 0, de + (j>)/dp + = 0, 
where p = S/Anr^ w 0.16 fm~ 3 is the density of the infinite nuclear matter and r = R/A 1 / 3 is the 
radius constant. Here we use a quadratic approximation e + = K(p + — p) 2 /(18p 2 ), where K is the 
incompressibility modulus of symmetric nuclear matter, mainly K m 220 — 260 MeV (see Table 
I). The isovector component can be simply evaluated as S- = J (I 2 - p 2 _/p\) [5]. The isoscalar 
SO gradient terms in Eq. (2) are defined with a constant: T> + = — 9mW / 2 /16/i 2 , where W ^100 - 
130 MeV-fm 5 and m is the nucleon mass [4, 23-26]. 

Minimizing the energy E under the constraints of the fixed particle number A = J dr p+(r) 
and neutron excess N — Z = J dr p_(r) (also others like a deformation [1, 3]) one arrives at the 
Lagrange equations with the corresponding multipliers, A + and A_ being the isoscalar and isovector 
chemical potentials, respectively (see Appendixes A and B). Our approach can be applied for any 
deformation parameter of the nuclear surface if its diffuseness with respect to the curvature radius 
a/R is small. The analytical solutions will be obtained approximately up to the order of A 2 / 3 in 
the binding energy. To satisfy the condition of the particle number conservation with the required 
accuracy we account for relatively small surface corrections (oc a/R ~ A^ 1 ^ 3 at the first order) to 
the leading terms in the Lagrange multipliers [2, 3, 5] (Appendixes B and C). 

For the isoscalar particle density, w = p + /p one has up to leading terms in the parameter a/R 
the usual first-order differential Lagrange equation with the solution [3, 5] 




(3) 



where /3 = T> + p/C + is the dimensionless SO parameter. For convenience we introduced also 
the dimensionless parameter e = 18e/K. For w r = w(x = 0) one has the boundary condition, 
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d 2 w{x)/dx 2 = at the ES (x = 0): 

e(w r ) + w r (l + /3uv) [de(w r )/dw] = . (4) 

In Eq. (3), a ^ 0.5 — 0.6 fm is the diffuseness parameter as shown in Table I [£ = (r — i?)/a for 
spherical nuclei in the spherical coordinates]. The diffuseness of the edge a d is given by 

a d = ^5(£2 -£ 2 )/3 = a^5(x~ 2 ~ -x 2 )/3 , x^ = I x n dx(dw/dx) , (5) 

J —CO 

where bars mean an averaging with the surface density distribution dw/dx [3]. For all Skyrme 
forces (Table I) the parameter a introduced in Eq. (3) measures the diffuseness of the nuclear 
edge as the mean-squared fluctuation of £ due to the relation a d a a/5/3 (see also Refs. [3, 15]). 
As shown in Ref. [5], the influence of the semiclassical h corrections (related to the ETF kinetic 
energy) to w(x) is negligibly small everywhere, besides the quantum tail outside the nucleus 
(x> 1). Therefore, all these corrections were neglected in Eq. (2). With a good convergence of 
the expansion of the e(y) in powers of 1 — y up to the quadratic term [3, 5], e = (1 — y) 2 , one 
finds the analytical solutions of Eq. (3) in terms of the algebraic, trigonometric and logarithmic 
functions [see Eq. (A3)]. For (3 = (i.e. without SO terms), it simplifies to the solution w(x) = 
tanh 2 [(x — xq)/2] for x < xq = 2arctanh(l/\/3) and zero for x outside the nucleus (x > xo). 

For the isovector density, w-(x) = p_/(pJ), after simple transformations of the isovec- 
tor Lagrange equation up to the leading term in a/R in the ES approximation one simi- 
larly finds the equation and boundary condition [see Eq. (A2)]. The analytical solution for 
u>_ = wcos[if)(w) / 'a/1 + P] can be obtained through the expansion (A5) of tp in powers of 

w(w) = (1 - w)/c sym with c sym = a ^ _ ^ | . (6) 

Expanding up to the second order in w one finds (Appendix A) 

w. = w(l- > = [1 + cw(w)} , c = ^7+ V 1 • (7) 

In Fig. 2 the SO dependence of the function w-(x) is compared with that of the density w(x) for 
the SLy7 force as a typical example [5]. In particular, as seen in Figs. 3 and 4 the isoscalar w(x) and 
the isovector w_ (x) densities depend only rather weakly on the most of the Skyrme forces [23, 24] . 
Only in a larger (logarithmic) scale, as shown in Fig. 4, one observes notable differences in the 
isovector densities W- derived from different Skyrme forces within the edge diffuseness. 



III. ISOVECTOR ENERGY AND STIFFNESS 



Within the improved ES approximation where also higher order corrections in the small param- 
eter a/R are taken into account we derive equations for the nuclear surface itself (see Appendix 
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B and Refs. [2, 3, 5]). For more exact isoscalar and isovector particle densities we account for the 
main terms in the next order of the parameter a/R in the Lagrange equations [cf. Eq. (Bl) as 
compared with Eq. (Al)]. Multiplying these equations by <9p_/<9£ and integrating them over the 
ES in the normal-to-surface direction £ and using the solutions for w±(x) up to the leading orders 
[see Eqs. (3) and (7)] one arrives at the ES equations in the form of the macroscopic boundary 
conditions (B2) [1-3, 5, 22, 28-30]. They ensure equilibrium through the equivalence of the volume 
and surface (capillary) pressure (isoscalar or isovector) variations. As shown in Appendix B, the 
latter ones are proportional to the corresponding surface tension coefficients: 

a ± = b^/^rl b ™ ~ MrlC ± jH d£ (l + ' . (8) 

The nuclear energy E (equation (1)) in this improved ES approximation (Appendix C) is split 
into volume and surface (both with the symmetry) terms, 

E w —b v A + J(N - Zf/A + E s . (9) 

For the surface energy Es one obtains 

E s = E s +) + E s - ] (10) 

with the following isoscalar (+) and isovector (-) surface components: 

E s ±) = a ± S = bf ) S/^rl) 1 (11) 

where S is the surface area of the ES. The energies E^ in Eq. (11) are determined by the isoscalar 
b s + ^ and isovector b s ~^ surface energy constants of Eq. (8). These constants are proportional to 
the corresponding surface tension coefficients o± in Eq. (8) through the solutions (3) and (7) for 
P±(0 which can be taken into account in leading order of a/R (Appendix C). These coefficients 
<t± are the same as found in the expressions for the capillary pressures of the boundary conditions 
(B2). 

For the energy surface coefficients b^ one obtains 



6 (+) = 6C +P j + ^ ^ = r 1 dwv / w(1 + /3w)e(w) ? (12) 

roa Jo 
b ( s ] = k s I 2 , k s = 6p C- J-/(r a) , (13) 



For w and c, see Eqs. (6) and (7), respectively. Simple expressions for the constants b s in Eqs. 
(12) and (13) can be easily derived in terms of algebraic and trigonometric functions by calculating 
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explicitly the integrals over w for the quadratic form of e(w) [Eqs. (C3) and (C4)]. Note that 
in these derivations we neglected curvature terms and, being of the same order, shell corrections. 
The isovector energy terms were obtained within the ES approximation with high accuracy up to 
the product of two small quantities, I 2 and (a/R) 2 . 

According to the theory [7-9], one may define the isovector stiffness Q with respect to the 
difference R n — R p between the neutron and proton radii as a collective variable, 

Ei^ = --fdSQr^-^ f , r = ^— , (15) 

where r is the neutron skin. Comparing this expression with Eq. (11) for the isovector surface 
energy written through the isovector surface energy constant [Eq. (13)], one obtains 

Q = -b^/r 2 = -k s I 2 /r 2 . (16) 

Defining the neutron and proton radii R n , p as the positions of the maxima of the neutron and 
proton density gradients, respectively, 

&2pnA -0, (^±) =0, (17) 



dr2 )^ Rn , v V dr 2 /r=R 



we use the expansion in small values of SR ntP = R n ,p — R near the ES. Thus, in the linear 
approximation in 5R UjP and I one obtains 



al d 2 w^ 



r dx 2 
where 



x=0 



d 3 w 



dx 3 



-l 



- {] J r oc 2 sym 



S " r g(w r ), (18) 



w 3 / 2 (l + Bw) 5/2 

9{W) = (l + (3)(3w + l + 4(3w) ( w ( 1 + 2 ^) 2 + 2^(l + c^)[c W -c s?/m (l + 2^)]}, (19) 

and w r is the solution of the boundary equation (4). In the derivations of Eq. (18), we used the 
approximation e(w) = (1 — w) 2 and expressions (3) for w(x) and (7) for w-(x)/w. The neutron 
and proton particle-density variations in Eq. (17) conserve the center of mass in the same linear 
approximation in 5R njP and /. Inserting Eqs. (18) and (13) into Eq. (16) one finally arrives at 

9 = -"k " = ^=«^)' (20) 

where J7_ and g(w r ) are given by Eqs. (14), (19) and (4). In the derivation of Eq. (20) we used 
also Eq. (3) for the diffuseness parameter a and Eq. (6) for the c sym . Note that Q = —9J 2 /4k s has 
been predicted in Refs. [7, 8] and therefore for v — 9/4 the first part of (20) which relates Q with 
the volume symmetry energy J and the isovector surface energy k s constants, is identical to that 
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used in Refs. [7-10, 14, 15]. However, in our derivations v deviates from 9/4 and it is proportional 
to the function J^ 2 /g 2 (w r ). This function depends significantly on the SO interaction parameter /3 
but not too much on the specific Skyrme forces. Indeed, the most sensitive parameter C_ cancels 
in the expression (20) for v: k s oc C_ and r oc l/c 2 sym oc C_ [see also Eqs. (13) for k s , (6) for c sym 
and (18) for r]. The constant v at /3 — can be easy evaluated, 

108 [l-8/(7c^)f 

25 l-4/(3c SJ/m ) ' V ; 

with neglecting small terms oc l/c 2 ym , c sym ~ 2 — 6 for the Skyrme parameters of Refs. [23, 24] 
[c sym = oo for T6 forces, see Eqs. (20), (14), (6) and (4) (w r = 1/3) and Table I]. Another 
difference in Q [ Eq. (20)] from that of Refs. [7-10] is the expression (13) itself for k s . Thus, 
the isovector stiffness coefficient Q introduced originally by Myers and Swiatecki [7, 8] is not a 
parameter of our approach but was found analytically in the explicit closed form (20) through the 
parameters of the Skyrme forces. 



IV. DISCUSSION OF THE RESULTS 



In Table II and also in Fig. 5 we show the isovector energy coefficient k s [Eq. (13)], the stiffness 
parameter Q [Eq. (20)] and the neutron skin r [Eq. (18)] obtained within the ES approximation 
using the quadratic approximation for e(w) for several Skyrme forces [23, 24] with parameters 
presented in Table I. We also show the quantities kso, v §i Qo an d t where the SO interaction is 
neglected (/3 = 0). One can see a fairly good agreement for the analytical isoscalar energy constant 
b^ (12) with that of Refs. [23, 24] and [5] (Table I). The isovector energy coefficient k s is more 
sensitive to the choice of the Skyrme forces than the isoscalar one (Eq. (12) and Ref. [5]). 
The modulus of k s is significantly larger for most of the Lyon Skyrme forces SLy [23] and SkI3 
[24] than for the other ones. For these forces the stiffnesses Q are correspondingly smaller. The 
isovector stiffness Q is even more sensitive to the constants of the Skyrme force than the constants 
k s . They are significantly larger for all forces, especially for SGII, than the well-known empirical 
values Q « 14 - 35 MeV [8-10]. 

Swiatecki and his collaborators [9] found the stiffness Q w 14 — 20 MeV by fitting the nuclear 
isovector giant dip ole- resonance (IVGDR) energies calculated in the simplest version of the hy- 
drodynamical model to the experimental data. Later, they suggested larger values Q ~ 30 — 35 
MeV accounting for a more detailed study of other phenomena in Refs. [8, 10]. In spite of several 
misprints in the derivations of the IVGDR energies in Ref. [9] (in particular, in Eq. (7.7) of [9] for 
the displacement of the center-of-mass conservation the factor NZ/A 2 should be in the numerator 
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and not in the denominator of the irrotational flow moment of inertia (see Ref. [31])), the final 
result for the IVGDR energy constant is almost the same as for the asymptotically large values of 
Q, ZJA'^/Q < 1 {NZ/A 2 « 1/4), 

D = hu_A 1/3 = ^oo/y / l + 3JA- 1 /3/g ! = sJ%h 2 J/(mrl) . (22) 

These values for D are in good agreement with the well-known experimental value D exp 80 
MeV for heavy nuclei (D 88 MeV) within a precision better or of the order of 20% 

(a little worse for the specific SkI3 Skyrme forces), as shown in Table II, see also Ref. [20] for a 
more proper HD approach and Refs. [32-34] for other semiclassical nuclear models taking all into 
account the nuclear surface motion. As shown in Ref. [6] the averaged IVGDR energies and the 
energy weighted sum rules (EWSR) obtained with the semiclassical FLD approach based on the 
Landau- Vlasov equation [22] with the macroscopic boundary conditions (see Appendix B) are also 
basically insensitive to the isovector surface energy constant ks, and they are similarly in good 
agreement with the experimental data. An investigation of the splitting of the IVGDR within 
this approach into the main peak which exhausts mainly the independent-of-model EWSR and 
a satellite (with a much smaller contribution to the EWSR) focusing on a much more sensitive 
ks dependence of the pygmy (IVGDR satellite) resonances (see Refs. [18, 19]) will be published 
elsewhere. 

More precise A-dependence of the quantity D [Eq. (22)] for finite values of Q seems to be 
beyond the accuracy of these HD calculations because of several other reasons. More realistic 
self-consistent HF calculations accounting for the Coulomb interaction, surface-curvature and 
quantum-shell effects led to larger Q pa 30 — 80 MeV [4, 15]. With larger Q (see Table II) the 
fundamental parameter of the LDM expansion in Ref. [7], (9J/AQ)A~ 1 ^ 3 , is really small for A> 40, 
and therefore, results obtained by using this expansion are more justified. 

The most responsible parameter of the Skyrme HF approach leading in significant differences 
in the ks and Q values is the constant C_ in the gradient terms of the energy density [Eq. (2) 
and Table I]. Indeed, the key quantity in the expression for Q, Eq. (20), and the isovector surface 
energy constant k s [or bg, Eq. (13)], is the constant C_ because one mainly has k s oc C_ (see Fig. 
5), and Q oc l/k s oc 1/C-. As seen in Table I and in Fig. 5 the constant C_ is very different for 
the different Skyrme forces (even in sign). As shown in Fig. 3 and below Eq. (20), other quantities 
in Eq. (20) are much less sensitive to most Skyrme interactions. The situation is very much in 
contrast to the isoscalar energy density constant b s + ^ oc C + [Eq. (12)]. All Skyrme parameters are 
fitted to the well-known experimental value b s + ^ = 17—19 MeV because C + is almost constant 
(Table I). Contrary to this, there are so far no clear experiments which would determine k s well 
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enough because the mean energies of the IVGDR (main peaks) do not depend very much on k s 
for the different Skyrme forces (see last two rows of Table II). Perhaps, the low-lying isovector 
collective states are more sensitive but there is no careful systematic study of their k s dependence 
at the present time. Another reason for so different ks and Q values might be traced back to the 
difficulties in deducing k s directly from the HF calculations due to the curvature and quantum 
effects in contrast to b^\ We have also to go far away from the nuclear stability line to subtract 
uniquely the coefficient k s in the dependence of oc I 2 = (N — Z) 2 /A 2 , according to Eq. (13). 
For exotic nuclei one has more problems to relate k s to the experimental data with a good enough 
precision. Note that k s is a more fundamental constant than the isovector stiffness Q due to 
the direct relation to the tension coefficient o"_ of the isovector capillary pressure. Therefore, it 
is simpler to analyze the experimental data for the IVGDR within the macroscopic HD or FLD 
models in terms of the constant k s . The quantity Q involves also the ES approximation for the 
description of the nuclear edge through the neutron skin r [see Eq. (15)]. The precision of this 
description depends more on the specific nuclear models [14-17]. On the other hand, the neutron 
skin thickness r is interesting in many aspects for the investigation of exotic nuclei, in particular, 
in nuclear astrophysics. 

We emphasize that for specific Skyrme forces there exists an abnormal behavior of the isovector 
surface constants k s and Q. It is related to the fundamental constant C_ of the energy density 
(2). For the parameter set T6 (C_ = 0) one finds k s = 0. Therefore, according to Eq. (20), the 
value of Q diverges [y is almost independent on C_). Notice that the isovector gradient terms 
which are important for the consistent derivations within the ES approach are also not included 
(C- = 0) into the symmetry energy density in Refs. [11, 13]. Moreover, for RATP [23] and Ski 
[24] (also for the specific Skyrme forces BSk6 and BSk8 of Ref. [12] 1 ), the isovector stiffness Q is 
even negative as C_ > (k s > 0) in contrast to all other Skyrme forces. 

Table II shows also the coefficients v of Eq. (20) for the isovector stiffness Q. They are mostly 
constant [u ~ 2 — 4, see Eq. (21) and Table II] for all Skyrme forces at ft — 0. However, these 
constants v are rather sensitive to the SO, i.e. to the /3 dependence of both the function g(w r ) 
[Eq. (19)] in the expression (18) for the neutron skin r and to the constant J- [Eq. (14)] in the 
expression (13) for the isovector energy coefficient k s . As compared to 9/4 suggested in Ref. [7], 
they are significantly smaller in magnitude for the most of the Skyrme forces (besides of SGII and 
T6 with larger values of v). 

1 Notice also that for the BSk forces [12] we found an unexpected behavior of the particle density W-(x) with a 
negative minimum near the ES (a proton instead of neutron skin because of p n < p p , in contrast to any other 
forces discussed in Ref. [23] with w- being always positive) 
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V. CONCLUSIONS 

Simple expressions for the isovector parts of the particle densities and energies in the leading ES 
approximation were used for the derivation of analytical expressions of the neutron skin thickness 
and the isovector stiffness coefficients. As shown in Appendix B we have to include higher order 
terms in the parameter a/R. These terms depend on the well-known parameters of the Skyrme 
forces. Results for the isovector surface energy constant ks, the neutron skin thickness r and the 
stiffness Q depend in a sensitive way on the choice of the parameters for the Skyrme functional, 
especially on the parameter C_ in the gradient terms of the density in the surface symmetry energy 
density of Eq. (2). The values of the isovector constants ks, r and Q depend also very much on the 
SO interaction constant f3. The isovector stiffness constants Q are significantly larger than those 
found earlier for all desired Skyrme forces. The mean IVGDR energies and sum rules calculated 
in the HD [9, 20, 31] and FLD [6, 22] models for the most k s values in Table II are in a fairly good 
agreement with the experimental data. For further perspectives, it would be worth to apply our 
results to the calculations of the pygmy resonances in the IVGDR strength within the FLD model 
[22] and the isovector low-lying collective states within the periodic orbit theory [35-37] which 
are expected to be more sensitive to the values of ks- Our approach is helpful for further study 
of the effects in the surface symmetry energy because it gives the analytical expressions for the 
constants k s and Q directly connected with a few critical parameters of the Skyrme interaction 
without using any fitting. 
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Appendix A: Solutions to the isovector Lagrange equation 

The Lagrange equation for the variations of the isovector particle density p_ in the energy 
density (2) up to the leading terms in a small parameter a/R is given by [5] 

C-%- + -£-\p + e-(p + ,P-)]=0, (Al) 
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where e_ is denned just below Eq. (2). We neglected here the higher order terms proportional 
to the first derivatives of the particle density p_ with respect to £ and the surface correction to 
the isovector chemical potential like in Refs. [2, 3] for the isoscalar case. For the dimensionless 
isovector density W- = p-/(pl) after simple transformations one finds the following equation and 
the boundary condition in the form 



dw_ 
dw 




w_(l) = 1 , (A2) 



where /3 is the SO parameter defined below Eq. (3), see also Eq. (6) for c sym . The above equation 
determines the isovector density W- as a function of the isoscalar one w(x) [Eq. (3)]. In the 
quadratic approximation for e(w) one explicitly finds 



x(w) = y/l + P In 



1 - w) (l + (1 + 2(3)w r + 2 V /(1 + /3)(1 + /M> r ) 



^l - w r ) (l + (1 + 20)w + 2y/(l + P)(l + l3w)iv) 
+ [arcsin (1 + 2/3w) - arcsin (1 + 2/3w r )\ , (A3) 

where w r is the solution of the boundary condition (4). Substituting w- = w cosip into Eq. (A2), 
and taking the approximation e = (1 — w) 2 , one has the following first order differential equation 
for a new function ip(w): 



w(l — w) . , ,, r — . , 1—w . 

sim^> ip — yl + pw smip cos^, ^(1J = 0. 



Csym Csym 



(A4) 



The boundary condition for this equation is related to that of Eq. (A2) for w-(w). This equation 
looks more complicated because of the trigonometric nonlinear terms. However, it allows to 
obtain simple approximate and rather exact analytical solutions within standard perturbation 
theory. Indeed, according to Eqs. (A2) and (3) where we have did not express the x-dependence 
explicitly we note that u>_ oc w(x) is a sharply decreasing function of x within a small diffuseness 
region of the order of one in dimensionless units (Figs. 2-4). Thus, we may find the approximate 
solutions to the equation (A4) (with its boundary condition) in terms of a power expansion of a 
new function ip(w) in terms of a new small argument w, 

oo 

$(w) = ip(w) = ^c n w n , (A5) 

?1=0 

with unknown coefficients c n and w defined in Eq. (6). Substituting the power series (A5) into 
Eq. (A4) one expands first the trigonometric functions into power series of w in accordance with 
the boundary condition in Eq. (A4). As usual, using standard perturbation theory, we obtain 
the system of algebraic equations for the coefficients c n [Eq. (A5)] by equating the coefficients 
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from both sides of Eq. (A4) at the same powers of w. This simple procedure leads to a system of 
algebraic recurrence relations which determine the coefficients c n as functions of the parameters 
(5 and c sym of Eq. (A4), 



c 3 



+ cl (c sym - ±) + l^TTpcl + c 2 - 3 Cl ) I , (A6) 



v/TT? [ 8 (l + /3) 3/2 

etc. In particular, up to the second order in w, we derive an analytical solution in an explicitly 
closed form: 

${w) = w (ci + c 2 w) , Cl = -j=L=, c 2 = + * (^c aym - lj . (A7) 

Thus, using the standard perturbation expansion method of solving tp(w) in terms of the power 
series of the w (up to w 2 ), one obtains the quadratic expansion of ip{w) [Eq. (7)] with c — c^/ci. 
Notice that one finds a good convergence of the power expansion of ip{w) (A7) in w for w-(x) at 
the second order in w because of the large value of c sym for all Skyrme forces presented in Table 
I [Eq. (6) for c sym \. 

Appendix B: The macroscopic boundary conditions and surface tension coefficients 

For the derivation of the expression for the surface tension coefficients cr±, we first write the 
system of the Lagrange equations by using variations of the energy density £(p+, p~) with respect 
to the isoscalar and isovector densities p + and p_. Then, we substitute the solution of the first 
Lagrange equation for the variations of the isoscalar density p = p + in the energy density (2) 
(Refs. [2, 3]) into the second Lagrange equation for the isovector density p_. Using the Laplacian 
in the variables £ and i] (Appendix A in Ref. [2]) we keep the major terms in this second equation 
within the improved precision in the small parameter a/R. The improved precision means that 
we take into account the next terms proportional to the first derivatives of the particle densities 
[along with the second ones of Eq. (Al)] and the small surface corrections A± to the isoscalar and 
isovector Lagrange multipliers \±. Within this improved precision, one finds the second Lagrange 
equation by the variations of the energy density £{p + ,p ), Eq. (2), with respect to the isovector 
particle density p_: 

C- ^ + K-H 9 -^- - ^- [p + e_( P+ , P-)\ + A_ = 0, (Bl) 

where H is the mean curvature of the ES (H = 1/R for the spherical ES). The isovector chemical- 
potential correction A_ was introduced [5] like the isoscalar one A + worked out in details in Refs. 
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[2, 3]. Multiplying Eq. (Bl) by <9p_/<9£ we integrate in the coordinate £ normal to the ES from a 
spatial point £j n inside the volume (at ^ — a ) to oo term by term. Using also the integration 
by parts, within the ES approximation it results in the macroscopic boundary conditions (together 
with the isoscalar condition from [2, 3, 5, 22, 29, 30]): 

(p I A-) ES = P]-) = 2a.H, (p A + ) ES = P« = 2a + H. (B2) 

Here, Pj 1 ^ are the isovector and isoscalar surface-tension (capillary) pressures and a± are the 
corresponding tension coefficients [see their expressions in Eq. (8)]. The Lagrange multipliers A± 
multiplied by pi and p in the circle brackets on the left hand sides of equations (B2) are the volume 
isovector (p/A_) and isoscalar (pA + ) pressure excesses, respectively (see Ref. [5]). These pressures 
due to the surface curvature can be derived using the volume solutions of the Lagrange equations 
for the particle densities [obtained by doing variations of the energy density S and neglecting all 
the derivatives of the particle densities in Eq. (2)], 

T QA , / 81 A , A 18. 7" J 

(B3) 



P- ~ P 



9A+\ A. 

1+ ir ) + 2j 





9A+ / 


, P+-P 


. 1+ K\ 



i _ 81A +A _ H£ p 

2K J 2K 



Inserting A + and A_ from Eq. (B2) into Eq. (B3) one gets 



p- = pi 



1 6b ( s +) H r 2b ( s ] H r 



(B4) 



K 6 J P 

As seen from Eq. (B4), the isovector density correction to the volume density p_ due to a finiteness 
of the coupled system of the two Lagrange equations depend on both isoscalar and isovector surface 
energy constants in the first order expansion of the small parameter a/R. If we are not too 
far from the valley of stability / is an additional small parameter and the isovector corrections 
are small as compared with the isoscalar values [b^ oc I 2 , A_ oc /, see Eqs. (B3), (B4) and 
(13)]. Thus, Eq. (B2) has a clear physical meaning as the macroscopic boundary conditions for 
equilibrium of the isovector and isoscalar forces (volume and surface pressures) acting on the ES 
[22, 28]. Note that in these derivations we used Eq. (Al) in leading order in a/R [neglecting small 
terms proportional to the first derivative <9p_/<9£ and A_ in Eq. (Bl)], like for the isoscalar case 
considered in Refs. [2, 3, 29, 30]. The obvious boundary conditions of disappearance of the particle 
densities and all their derivatives with respect to £ outside of the ES for £ — > oo (£ 3> a) were 
taken into account too. The isovector tension coefficient cr_ is much smaller than the isoscalar 
one cx_|_ [see Eq. (8)] as ct_ oc I 2 due to p_ oc / and I <^ 1 is small near the nuclear stability line. 
Another reason is also because of a smallness of C_ as compared to C + for the realistic Skyrme 
forces [23, 24]. ^From comparison of Eqs. (B3) and (B4) for p_ [see also Eq. (8)], one may evaluate 

A o_ rr ol(-) 
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which is consistent with Eq. (Bl) (r H ~ a/R in these estimations, see corresponding ones in 



For calculations of the surface energy components Eg of the energy E, Eq. (1), within the 
same improved ES approximation as described above in Appendix B we first may separate the 
volume terms related to the first two terms of Eq. (2) for the energy density £ . Other terms of 
the energy density £(p + ,p_) in Eq. (2) lead to the surface components E^ [Eq. (11)], as they 
are concentrated near the ES. Integrating the energy density £ [see Eq. (2)] over the spatial 
coordinates r in the local coordinate system £,77 (see Fig. 1) in the ES approximation one finds 



where — a is like in Appendix B [2, 3, 5]. The local coordinates £,r? were used because 

the integral over £ converges rapidly within the ES layer which is effectively taken for [£[< a. 
Therefore, we may extend formally £ iri to —00 in the first (internal) integral taken over the ES 
in the normal direction £ in Eq. (CI). Then, the second integration is performed over the closed 
surface of the ES. The integrand over £ contains terms of the order of (p/a) 2 oc (R/a) 2 like the 
ones of the leading order in Eq. (Bl) [see for instance the second derivatives in Eq. (Al) which 
are also oc (R/a) 2 }. However, the integration is effectively performed over the edge region of the 
order of a that leads to the additional smallness proportional to a/R. At this leading order the 
rj dependence of the internal integrand can be neglected. Moreover, from the Lagrange equations 
at this order one can realize that the terms without the particle density gradients in Eq. (CI) are 
equivalent to the gradient terms. Therefore, for the calculation of the internal integral we may 
approximately reduce the integrand over £ to the only derivatives of the particle densities of the 
leading order p±(£) in £ (with the factor 2) using (Vp±) 2 — p+e± (p + ,p_) ~ 2(<9p±/<9£) 2 [see Eqs. 
(3) and (7)] . Taking the integral over £ within the infinite integration region (—00 < £ < 00) off 
the integral over the ES (d«S) we are left with the integral over the ES itself that is the surface 
area S. Thus, we arrive finally at the r.h.s. of Eq. (CI) with the surface tension coefficient a± 



Using now the quadratic approximation e(w) = (1 — w) 2 in Eq. (8) for b$ = AKr 2 .a± (V_ = 0) 
one obtains (for f3 < 0, see Table I) 



Refs. [2, 3]). 



Appendix C: Derivations of the surface energy and its coefficients 




(CI) 



[Eq. (8)]. 




6p C± J±/(r a) , 



(C2) 
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where 



J, 



= / dw \/vj( 
Jo 



(2) 

+ J\ arcsin a/ — p 



jj 1 ) = 3 + 4/3(1 + /?), j7^ = -3-6/3 



(2) 



For the isovector energy constant J_ one finds 
1 Z" 1 

X = -- 3 / dw ^w{l + Pw) (l-w)(l + cw) 2 

1 + P Jo 

J { - ] {csyjc) y/-0(l + /3) + J {2) (c syrn /Z) arcsin 



c 2 



1920(1 + /3)(-/3) 9 / 2 



^(C) = 105 - 4/3 {95 + 75C + /3 [119 + 10C(19 + 6() + 8/3 2 (1 + 10C(1 + 0) 



+ 8C(5C(3 + 2C)-6)]}, 



j[ 2 \() = 15 {7 + 2/3 [5(3 + 2C) + 8/3(1 + C) (3 + C + 2/3(1 + C))]} . 



(C3) 



-0 



(C4) 



These equations determine explicitly the analytical expressions for the isoscalar (bg) and isovector 
(bg) energy constants in terms of the Skyrme force parameters, see Eqs. (7) for c, (6) for c sym and 
w. For the limit /3 ->■ from Eqs. (C3) and (C4) one has J± ->■ 4/15. With Eqs. (18) and (19) 
one arrives also at the explicit analytical expression for the isovector stiffness Q as a function of 
C- and /3. In the limit C_ — > one obtains fc s — )■ and Q — > oo because of the finite limit of the 
argument c sym /c — )• (1 + /3)//3 of the function J71 in Eq. (C4) [see also Eqs. (7) for c and (6) for 
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FIG. 1. The ES and local £, r/ coordinates. The ES in the cylindrical y, z coordinates with the symmetry 
axis z and diffuseness parameter a are shown schematically by the thick solid and dashed curves (after 
Ref. [5]). 
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SkM* 


SkM 


SIII 


SGII 


RATP 


SkP 


T6 


SkI3 


SLyb 


SLy4 


SLy5 


SLy6 


SLy7 


p, fm 6 


0.16 


0.16 


0.15 


0.16 


0.16 


0.16 


0.16 


0.16 


0.16 


0.16 


0.16 


0.17 


0.16 


b v , MeV 


15.8 


15.8 


15.9 


15.6 


16.0 


15.9 


16.0 


16.0 


16.0 


16.0 


16.0 


16.0 


16.0 


K, MeV 


217 


217 


355 


215 


240 


201 


236 


258 


230 


230 


230 


230 


230 


J, MeV 


30.1 


31.0 


28.2 


26.9 


29.3 


30.0 


30.0 


34.8 


32.0 


32.0 


32.1 


31.9 


31.9 


C+, MeVfm 5 


57.6 


52.9 


49.4 


43.9 


60.2 


60.1 


55.1 


51.8 


59.5 


59.5 


59.3 


54.0 


52.7 


C-, MeVfm 5 


-4.79 


-4.69 


-5.59 


-0.94 


13.9 


-20.2 





12.6 


-22.3 


-22.3 


-22.8 


-15.6 


-13.4 


a fm 


0.52 


0.50 


0.59 


0.45 


0.55 


0.50 


0.52 


0.53 


0.53 


0.53 


0.53 


0.52 


0.50 


a d fm 


0.63 


0.58 


0.73 


0.58 


0.71 


0.71 


0.70 


0.63 


0.68 


0.68 


0.68 


0.61 


0.61 


C.sym 


3.26 


3.21 


3.42 


6.02 


2.00 


1.52 


oo 


2.20 


1.59 


1.59 


1.57 


1.80 


1.93 


P 


-0.64 


-0.69 


-0.57 


-0.54 


-0.52 


-0.37 


-0.45 


-0.65 


-0.55 


-0.55 


-0.58 


-0.59 


-0.65 


b%, MeV [23, 24] 


16.0 


16.0 


17.0 


14.8 


17.9 


17.9 


17.9 


16.0 


16.7 


18.1 


18.0 


17.4 


17.0 


6+, MeV 


21.2 


19.9 


14.5 


18.7 


21.7 


24.9 


21.3 


18.3 


21.7 


21.7 


21.5 


21.6 


19.6 



TABLE I. Basic parameters of the Skyrme forces from Refs. [23, 24] and the isoscalar surface energy 
constants of Eq. (12); the critical parameters a [Eq. (3)] and ad [Eq. 5] for the nuclear diffuseness 
edge, the isoscalar and isovector constants C± of the energy density [Eq. (2)], c sym [Eq. (6)] and the 
spin-orbit constant (3 [see below Eq. (3)]; SLyb denotes shortly SLy230b of Ref. [23]. 
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FIG. 3. Isovector particle densities w-(x) (7) as functions of x within the quadratic approximation to 
e(w) for several Skyrme forces [23], see also Ref. [5]. 
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FIG. 4. The same as in Fig. 3 for the critical Skyrme interactions but within a smaller edge diffuseness 
region in the logarithmic scale. 
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FIG. 5. Isovector energy constant k s Eq. (13) (solid bars, left axis) vs the coefficient C_ (dotted points, 
right scale) of Eq. (2); k s fl (dashed bars) is given by Eq. (13) without the spin-orbit interaction (/3 = 0). 



22 





SkM* 


SkM 


SIII 


SGII 


RATP 


SkP 


T6 


SkI3 


SLyb 


SLy4 


SLy5 


SLy6 


SLy7 


k s ,o, MeV 


-3.26 


-3.84 


-2.65 


-0.71 


5.25 


-5.36 





6.93 


-7.51 


-7.54 


-8.14 


-7.45 


-6.95 


k s , MeV 


-0.77 


-1.90 


-0.52 


-0.21 


1.42 


-1.93 





4.88 


-4.24 


-4.38 


-6.96 


-6.72 


-6.32 


^0 


3.08 


3.06 


3.14 


3.64 


2.38 


2.17 


4.32 


2.53 


2.12 


2.12 


2.12 


2.22 


3.32 


V 


0.34 


0.46 


1.42 


17.9 


0.45 


1.76 


4.30 


0.56 


0.44 


0.44 


0.59 


0.65 


0.67 


Qo, MeV 


7744 


9487 


6255 


16879 


-371 


3815 


oo 


-6314 


4771 


4794 


5178 


5350 


5703 


Q, MeV 


398 


234 


2168 


60998 


-270 


823 


oo 


-140 


105 


104 


87 


98 


109 


To// 


0.021 


0.020 


0.021 


0.0065 


0.038 


0.037 





0.033 


0.040 


0.40 


0.040 


0.037 


0.035 


r/J 


0.044 


0.090 


0.015 


0.0019 


0.072 


0.048 





0.187 


0.20 


0.21 


0.28 


0.26 


0.24 


Dhd, MeV 


85-86 


85-86 


82 


82 


90-89 


87 


88 


105-100 


81-84 


81-84 


79-83 


81-85 


81-84 


Dfld-, MeV 


73-82 


71-76 


79-104 


74-77 


77-87 


70-69 


86-88 


101-106 


80-90 


80-90 


76-84 


80-91 


77-89 



TABLE II. Isovector energy k s and stiffness Q coefficients are shown for several Skyrme forces [23, 24]; 
v is the constant of Eq. (20); t/I is the neutron skin calculated by Eq. (18); quantities ks : o, v 0, Qo an d 
To are calculated with (3 = 0; the intervals of monotonic functions Dhd(A) and Dfld(A) for the HD 
and FLD models in the last two lines are related to A « 50 — 200 (the last line is taken from Ref . [6] ) . 



